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Abstract 

A perfect fluid is quantized by the canonical method. The con- 
straints are found and this allows the Dirac brackets to be calculated. 
Replacing the Dirac brackets with quantum commutators formally 
quantizes the system. There is a momentum operator in the denom- 
inator of some coordinate quantum commutators. It is shown that it 
is possible to multiply throughout by this momentum operator. Fac- 
tor ordering differences can result in a viscosity term. The resulting 
quantum commutator algebra is 



W4(«3«2 - V2V3) = -i, 

w 4 (uiu 3 - V3V1) = -iv 3 , 
V4V1 — V1V4 = —i, 
V5V2 - v 2 v 5 = -i. 
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1 Introduction 



It has been known for some time Hargreaves (1908) |IJ that a perfect fluid 
has a Lagrangian formulation. The Lagrangian is taken to be the pressure 
and variations are achieved through an infinitesimal form of the first law of 
thermodynamics. A perfect fluid's stress is described using the vector field 
co-moving with the fluid; this vector field defines an absolute time for the 
system, furthermore this absolute time can then be used to define canonical 
momenta and canonical Hamiltonians. This is done here for the first time. 
There are equivalences between scalar fields and fluids, Tabensky and Taub 
(1973) more generally the co-moving vector field can be decomposed 
into scalar fields resulting in a description of a perfect fluid employing only 
scalar fields. Previously this decomposition has been investigated by choos- 
ing an ad hoc global time rather than absolute time and defining canonical 
momenta and other quantities with respect to the global time. Typically the 
resulting theory is applied to cosmology Tipler (1986) ||] Lapchiniskii and 
Rubakov (1977) 0]. Once the constrained Hamiltonian has been calculated 
by the standard canonical method Dirac (1963) || (see also Hanson, Regge, 
and Teitelboim (1976) fl(|), the Dirac bracket can be replaced by quantum 
commutators. The original motive for investigating this was to find a fluid 
generalization of the Higg's model Roberts (1997) IJ. A quantum treatment 
is required to estimate the VEVs (quantum vacuum expectation values) of 
the scalar fields, the VEVs are related to the induced non-zero mass. The 
quantum commutator algebra is unusual, perhaps reflecting the structure of 
the scalar field decomposition of the co- moving vector field Eckart (1960) 
§], Selinger and Whitham(1968) §, and Rund (1979) g. It is hoped that 
eventually the present theory will be applied to low temperature super flu- 
ids. To do this it probably will be necessary to include a chemical potential 
term in the first law of thermodynamics as expressed by equation p]. 
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2 Lagrangian and Hamiltonian Formulation of a 
Perfect Fluid's Dynamics. 



A perfect fluid has a Lagrangian formulation in which the Lagrangian is the 
pressure p. Variation is achieved by using the first law of thermodynamics 

dp = ndh — nTds, (1) 

where n is the particle number, T the temperature, s the entropy, and h the 
enthalpy. The pressure and the density are equated to the enthalpy and the 
particle number by 

p + jjL = nh. (2) 

In four dimensions a vector can be decomposed into four scalars, however 
the five scalar decomposition 

a=2 

hVa = W a = <Ta+Y, 9 " a «' V°V a = -1, (3) 
a=l 

where a,b,c . . . are spacetime indices and a, j3, 7 ... are fluid scalar potential 
indices. For a = 1, s and 8 = j Tdr have interpretation as the entropy and 
the thermasy respectively, for a = 2 there is no such interpretation and no 
"second tempreture" as #2 = 0. From now on the index a is suppressed 
as it is straightforward to reinstate. There are other conventions for this 
scalar field decomposition, for example with a minus(-) instead of a plus(+) 
before the summed fields, "q" is used to notate an arbitrary scalar field, 
i.e. q = a,6, ors. The coordinate space action is taken to be 

I = J \f~9 P dx 4 . (4) 

Replacing the first law with dp = —nV a dW a — nTds, variation with respect 
to the metric gives 

Tab = (p + nWaYb + pg a b, (5) 
variation with respect to a, 6, and s gives 

{nV a ) a = h + nQ = 0, s = 0, 6 = T, (6) 

respectively. = V® is the expansion of the vector field. 

The canonical momenta are are given by IT 1 = ^ and arc 

IF = -n, U e = 0, IT = -n9. (7) 
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The Hamiltonian density is usually defined in terms of components of the 
canonical stress as Q\. In the present case the canonical stress is not defined 
so that the metric stress TP 1 is used instead; also 4-vectors are used rather 
than components, resulting in 

H d = V a V b T ab = (i. (8) 

The standard Poisson bracket is 

r , „ 5 A SB SA SB 

where i, which labels each field, is summed; i differs from a in that it sums 
over all fields for example a, 8\, s\, 62, S2 in four dimensions, a labels the 
sets 0i, si and 02, S2- The integral sign and measure have been suppressed 
and the variations are performed independently. When absolute time is used 
Hamiltons equations have an additional term in the expansion []l3| explicitly 



where H c is the canonical Hamiltonian 

H c = [ H d ^dx\ (11) 



From [7| the momenta are constrained 

<Pi=U s - 8U a , 2 = IT 9 . (12) 
The initial Hamiltonian is 

Hi = Utq - C, (13) 

replacing the dependent Tl's gives the ordinary hamiltonian 

H = W (<j + 0s) - C, (14) 

adding the constraints gives the constrained Hamiltonian 

H x = H + \ a <j) a , X l = s, X 2 = 8, (15) 

where the X's are the Lagrange multipliers. Replacing the ordinary Hamil- 
tonian in this gives the Hamiltonian density 

H d = W (a + 0s) + A 1 (n s + 0W) + A 2 n e - C, (16) 
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Substituting the values of the momentum the Hamiltonian density is still 
weakly the fluid density. Using Hamiton's equations in the form 10 the time 
evolution of any variable X is given by 

replacing the Hamiltonian density by H\ and then holding the multipliers 
constant so that 

{X,H X } = {X,H } + \ L {X,4> L }, (18) 
where l,k, . . .ranges over the multipliers, gives the time evolution 

dr or oil 



dX /. .i,\SX xl SX ~5X 

8X 



+ ((V a W 7 ) a — GH 

+ ((A 1 - s)ir - en 
+ ({y a eu a ) a - err) 



5X 



5n e 

SX 



STL 8 



ax . 5X -SX . <5X ,„ .. 5X 

" * + ^ + V-"sf- (, '" ) sf- (19) 

where ~ means is 'weakly equal to'. Letting X equal the constraints gives 
^r 1 ~ 0, this shows that there are no further constraints so that the Dirac 

dr 

brackets can now be constructed. 
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A quantity R(q,H) is first class [10] if 

{J2,&}ra0, 1 = 0,1, 



(20) 



otherwise it is second class. The second class constraints encapture the way 
in which there are more variables describing the sytem than is necessary, 
they give rise to the Dirac matrix C iR , c.f. ju| page 10, defined by 



C UK = {4>t,,(j) K }. 
In the present case it is given by 

' - 1 \ 



(21) 



c. „ 



IF x 



1 

IF 



1 J 

' - 1 
1 



-ia 2 U a , 



+i 



a 



if' 



(22) 



where C Lt } is the inverse of C iK and a 2 is a Pauli matrix, Bjorken and Drell 
(1965) flu) page 378, 



- i 
+i 



The Dirac bracket is defined by c.f. [1C] page 11 

{A, By* = {A, By - {A, fayc-Hfa, B}. 

In the present case this gives the Dirac bracket 

iA, B y* = iA, B y + ± 5 4(^-e s -A + n*m 



(23) 



(24) 



IF 



5s 5a 



n CT 56 \5s 5a 5Ii e 



Consistency is checked by noting 

H X = H -iH ,^yCZ K \ 

x K = -{H ,pyc-\ 

from which H\ given by ^ can be recovered with the correct A 's. 



(25) 



(26) 
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3 Quantization. 



To quantize a classical dynamical system the Dirac bracket is replaced by 
the commutator 



1 r 



{A, By* 

in 



AB - BA 



(27) 



where Ti is Planck's constant divided by 2ir the hat on A and B signifies 
that the variables are now operators, from now on the hat is dropped as 
the ellipis is immediate. There are various correspondence criteria which 
one would hope to investigate, for example: as fi — > there should be 
firstly the same time evolution, secondly the same stress, thirdly the first 
law should be recovered. Another correspondence criteria can be called the 
particle number criteria: the particle number n (see equation ||) should bear 
a relation to the quantum particle number constructed from creation and 
destruction operators. An intermediate aim, between formal quantization 
achieved by replacing field and momenta Dirac brackets with commutators, 
and establishing contact with applications is to produce a quantum perfect 
fluid. This could be obtained from brackets involving the numbered field, 
the angular momentum and so on, or from brackets involving a mixture 
of these and geometric objects. However no progress has been made so 
far in finding a quantum perfect fluid, so that attention is restricted to the 
implications of replacing brackets consisting solely of individual components 
of fields and momenta with commutators. Effecting the replacement of the 
15 Dirac brackets between the fields and momenta there are 4 non-vanishing 
commutators 

[ Wo - o-w] = -m, [we - en 9 ] = o, [Ws - sW = -in, 

9 1 

[ a e -9o-} = -ih—, [9s-s9] = -ih—, [as-sa]=0. (28) 

The last two commutators have the operator W in the denominator. This 
might not be well-defined. To avoid W in the denominator we multiply 
by the operator W, using the first commutation of ^ it turns out that 
multiplying on the left or multiplying on the right are equivalent so that 

[aO - 9a] W = W [ad - 9a] = -i%9, 

[9s - s9] W = W [9s - s9] = -ih. (29) 

These results are in accord with the equations deduced if the Dirac brack- 
ets {.q 1 ,qill k y* are replaced by commutators. Left and right multiplication 
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by n CT are also equivalent if anti-commutation rather than commutation is 
considered. 

The quantum Hamiltonian is 



H q = + hU a a + l 2 aIF 
- l 5 Ii a s-l G 9sIl 



l 3 U a 9s - Z 4 lTs0 
- l 7 sW6 - l 8 s9U e 



p 



(30) 



where the l's are constant and obey h + h = 1 and h + h + h + le + h + ls = !• 



using the commutation relations the quantum Hamiltonian 30 becomes 



ihOl 



(31) 



where I = I2 + h + h + is called the ordering constant: it is of undefined 
size but is can be taken to be of order unity. Because the Dirac bracket of 
p with anything vanishes the commutators with p also vanish and p can be 
taken to be pi, where 1 is the identity element. To investigate the algebraic 



implications of 2£ and 31 label the 6 operators by v's, 

Vl = a,v 2 = s, v 3 = e,v 4 = rr, v 5 = u s ,v 6 = n e . (32) 

vq commutes with everything and can be disregarded. Of the remaining 



commutators only 4 are non-zero. In units H = 1, 29 and 32 give the algebra 

Vi(viv 3 - v 3 vi) = -iv 3 , 



V4V1 — V1V4 
V5V2 — V2V5 



(33) 



This algebra does not seem to be realizable in terms of matrices and differ- 
ential operators, the closest algebras are found in Jl2|| . If a commutator is 
constructed with a time derivative of the field or momenta, the same alge- 
bra results but multiplied by a term in the expansion. Similarly if m time 
derivatives occur, the algebra is multiplied by the expansion to the power 
of m. 
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